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ABSTRACT 

It has been shown that F, G, and early K dwarf hosts of Neptune-sized planets are not preferentially 
metal-rich. However, it is less clear whether the same holds for late K and M dwarf planet hosts. We 
report metallicities of Kepler targets and candidate transiting planet hosts with effective temperatures 
below 4500 K. We use new metallicity calibrations to determine [Fe/H] from visible and near- infrared 
spectra. We find that the metallicity distribution of late K and M dwarfs monitored by Kepler is 
consistent with that of the solar neighborhood. Further, we show that hosts of Earth- to Neptune- 
sized planets have metallicities consistent with those lacking detected planets and rule out a previously 
claimed 0.2 dex offset between the two distributions at 6cr confidence. We also demonstrate that the 
metallicities of late K and M dwarfs hosting multiple detected planets are consistent with those lacking 
detected planets. Our results indicate that multiple terrestrial and Neptune-sized planets can form 
around late K and M dwarfs with metallicities as low as 0.25 of the solar value. The presence of 
Neptune-sized planets orbiting such low-metallicity M dwarfs suggests that accreting planets collect 
most or all of the solids from the disk and that the potential cores of giant planets can readily form 
around M dwarfs. The paucity of giant planets around M dwarfs compared to solar-type stars must 
be due to relatively rapid disk evaporation or a slower rate of core accretion, rather than insufficient 
solids to form a core. 

Subject headings: Planetary Systems — planets and satellites: detection — planets and satellites: 
formation — stars: late-type — stars: abundances 



1. INTRODUCTION 



The NASA Kepler mission (jBorucki et al.ll2010( ) has 
discovered more than 2000 exoplanet candidate s (also 
called Kepler Objects of Interest or KOIs, Batalh a et al.l 
|2013() . enabling the study of exoplanet statistics based 
on large data sets. Among other science results, Ke- 
pler data has been used to estimate planet occurrence 
(e.g.. iHoward et al.ll2012t iFressin et a l."2013V constrain 
the distribution of planet d ensities {Gaidos et al. 2012; 
iWolfgang fc LaughlinI 20121). study the a rchitecture of 



multi-planet systems ( Fabrvckv et al.|[2012;il . and search 
for correlations (or non-correlations) between the radius 
of planets and the met allicity of their host s tars for F, 
G, and early K dwarfs (jBuchhave et al.|[2012[) . 

It is well established that the presence of Jovian plan- 
ets is correlated with the me t allicity of the host star 
(iGonzalezI 119971: iSantos et all 120011: iFischer fc Valentil 
I2005D . This is generally interpreted as supporting the 
core accretion mechanism for giant planet formation; 
metal-rich stars are assumed to have had metal-rich disks 
in which the higher density of solids allowed faster core 
accretion and the formation of gia n t plan ets before the 
disks dissipated. iBuchhave et al.l (|2012D showed that 
Earth- to Neptune-sized planets are present around FGK 
dwarfs spanning a ran ge of metallic i ties (— 0.6 < [Fe/H] < 
-fO.5). However, Bu chhave et al.l (|2012[ ) did not mea- 
sure the metallicity of a control sample of Kepler field 
stars. If Kepler is biased towards more metal-poor stars 
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(compared to the solar neighborhood), t he stars hosting 
Neptu ne-sized planets investigated by IBuchhave et al.l 
(|2012f ) c ould be more metal-ri ch than non-hosts. Fur- 
ther, the IBuchhave et aTl ()2012n sample contains no stars 
with Toff < 4500 K, and cannot draw conclusions about 
the role of metallicity on the frequency of Neptune-sized 

p lanets around la t e K a nd M dwarfs. 

iLaughlin et"aI1 (|2004[ ) and lAdams et aH (|2005f ) ar- 
gued that the core-accretion model of planet forma- 
tion predicts that late K and M dwarfs have sig- 
nificantly fewer giant planets than their solar-type 
counterparts. Disks around M dwarfs have longer 
dynamic al (orbital) times ( resulting in slower planet 
growth, iPollack et al.l 119961 ). lower surface densities 



(jHartmann et al.lll998l: LS cholz e t al.1 120061) . and less to- 
tal disk mass (jWilliams fc Cicza 201 1[) than those ar ound 
their solar-m ass counterparts. iLaughlin et al.l ()2004[ ) and 
lAdams et a l. (2005) predicted that although M dwarfs 
should have fewer giant planets, Neptune-like objects and 
terrestrial-type planets should be common around such 
stars. 

There are theoretical reasons to suspect the presence 
of Ncptunes around M dwarfs should be correlated with 
stellar metallicity, even if this is not the case for FGK 
dwarfs. Numerical simulations indicate that the ini- 
tial surface density of solids in a disk (for which stellar 
metallicit y is a proxy) c o ntrols the mass and number of 
planets. iKokubo et al.l (|2006f ) found that the mass of 
the largest and second largest planet in a planetary sys- 
tem should scale almost linearly with the disk surface 
density and that the total number of planets decreases 
with surface density, even in the absence of giant plan- 
ets. Because the surface density of solids in a planet- 
forming disk should scale with the metallicity, their re- 
sults suggest that metal-rich systems host larger (non- 
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Jovian) but fewer planets. Since the disk mass scales 
roughly linearly with the stellar mass (alth ough with con- 
siderable scatter: [Williams fc Ciezall201lD . it is possible 
that even metal-poor FGK dwarf disks have sufficient 
solid material to produce Neptune sized planets, as was 
found qbservationallv ( Sousa et al. 2008: Buchhave et alJ 
|2012[) . iSchlaufman fc Laughhnr()2010[ ) asserted that be- 
cause of the smaller disk mass around late K and M 
dwarfs, metallicity is more critical for the formation of 
Neptunes around these stars. They claimed this require- 
ment should manifest itself as a correlation between the 
presence of Ncptune-sized planets and the metallicity of 
late K or M dwarfs. 

Indeed, ISchlaufman fc LaughlinI ()201lD found that late 
K and M dwarfs hosting super-Earth- to Neptune- 
sized transiting planet candidates from Kepler have red- 
der g — r (for a fixed J — H) color than those with 
no detected transit. Based on a comparison between 
two open stellar clusters wit h different metallicities, 
ISchlaufman fc LaughlinI (|2011l ) claimed that the g — r 
color offset is due to a difference in metallicity of ~ 
0.2 dex, in agreement with the t heoretical case laid out 
in Schlaufman & Laug hlml (|2010l ). 

Complicating t he i s sue, the clusters used by 
ISchlaufman fc LaughlinI (|2011f) to calibrate their color- 
metallicity r elation cont ain very few late K and M stars. 
iWest et al.l (|2004D and iBochanski et al.l (|2013D found 
that the metallicity dependence oi g — r reverses sign at 
late K/early M spectral types, and that cool er stars have 
bluer a — r colors if they are more metal rich. iMann et al.l 
(|2012l ) explained that the origin of the g — r color differ- 
ence between the KOI and non -KOI sample observed by 
ISchlaufman fc LaughlinI ()2011[ ) is an artifact of giant star 
contamination in their non-KOI sample. But the ques- 
tion of whether KOI M dwarfs are more metal-rich than 
non-KOIs remains open. 

Compared to solar-type stars, M dwarf metallicities 
are difficult to determine, primarily due to the presence 
of complex molecular lines in their visible spectra, which 
result in line confusion and a lack of identifiable contin- 
uum, a nd do not always m atch with current M dwarf 
models (jAUard et al.ll20li[) . Previous techniques to de- 
termine M dwarf metallicities using color-magnitude di- 
agram s (| Johnson fc AppsI 120091 : ISchlaufman fc LaughlinI 
|2010| ) require astrometric parallaxes, which are largely 
unavailable for Kepler targets. Visible-light spectro- 
scopic techniques were developed ( e.g., the C index, 
iLepine et al]l2007l : JDhital et al.ll2012[ ). and are a reliable 
indicator of whether an M dwarf is a sub-dwarf or ul- 
tra su b-dwarf. However , C saturates near solar metal- 
hcity (jWoolf et al.' '2009': 'Mann et al.l l2013D . making it 
unreliable f or m easuring metallicities hig her than that 
of the Sun. iRojas-Avala et al.l ()2010l ) and iTerrien et al.l 
(|2012f ) demonstrated that atomic lines in the K- and 
iJ-bands (respectively) can be used to estimate metal- 
licities for M dwarfs. They calibrated their methods 
using ^20 wide binaries, but their samples were re- 
stricted in both spectral type (M0- M4), and metallicity 
(-0.5 < [Fe/H] < -hO.4). Recently, IMann et al.l (poTl 
(henceforth M13) used 110 wide binaries spanning K5 to 
M6, and —1.04 < [Fe/H] < -1-0.56, and derive improved 
calibrations to determine metallicities using visible, J-, 
H-, or i^-band spectra. 

Using the techniques of iRoias-Avala et afl (|2012D . 



iMuirhead et al.l ()2012af ) estimated met allicities for late 
K and M planet candidate hosts. IMuirhead et al.l 
(|2012aD found that late K and M KOIs' metallici- 
ties ar e consistent with the solar neighborhood (~ 
—0.10, iCasagrande et al.l 120081 ). but did not measure 
the metallicity of the o verall field for comparison. 
iDressing fc Charbonneaul (|2013f ) fitt ed qrizJHK colors 
from the Kepler In put Catalog (KIC. lBatal ha et al. l2010l; 
iBrown et al.l 120 lU) to the stellar models of Potter et al.l 
(J200S) to determine K,,M^, Teff, and [Fe/H] for M dwarf 
Kepler targets. However, the colors of late K and early- 
M dwarfs are usually not reliable indic ators of metallicity 
(ILepine et al.l 120131: IMann et al.l 2013 ) , and metallicities 
from IDressing fc Charbonnead ( 2013 ) s how weak or no 
significant correlation with those from IMuirhead et al.l 
(l2_Qi2a!). 

In this paper we investigate wheather the size and mul- 
tiplicity of planets around late-type dwarfs depend on 
the metallicity of the host star. In Section |2| we describe 
our sample of planet candidate hosts and our comparison 
sample of dwarfs with no detected transit. In Section [3] 
we detail our visible and near-infrared observations of 
Kepler stars. We derive a new calibration in Section |4| 
to determine [Fe/H] from visible wavelength spectra. We 
then apply this calibration, and others from M13, to cal- 
culate the [Fe/H] of the KOI and non-KOI samples. In 
Section[5]we report the metallicity distributions of late K 
and M dwarf hosts of Earth-, Neptune-, and Jovian-sized 
planets, hosts of multiple detected planets, and dwarfs 
with no detected transits. In Section |6| we conclude with 
a brief discussion of possible complications and the con- 
sequences of our findings. 

2. SAMPLE 

2.1. KOI Sample 

We selected KOIs from iBatalha et all (|2013D with 
Kp — J > 1.85, where Kp is the magnitude in the Ke- 
pler bandpass, a nd J is from the Tw o Micron All Sky 
Survey (2MASS. ISkrutskie et aLll2006( ). This sample in- 
cludes all dw arfs with T^ff < 4100 K and some as warm 
as 4500 K (Ma nn et al.1l2012D . We excluded 4 KOIs that 
are probably false posit ives. KOI-977 is a a giant star 
(IMuirhead et al.ll2012al) . KOL1902 has a V-like tran- 
sit shape, and flux variations indicative of an eclipsing 
binary. In the latest planet candidate release, KOI-1164 
has been added to the false positive lisO The light curve 
of KOI-256 shows no limb-darkening, which is indicative 
of a stellar eclipse, rather than a planetary transit. As a 
test, we obtained two spectra of KOI-256 6h apart (see 
Section |3| for a description of observations) . The spec- 
tra show a radial velocity difference of > 100 kms^^, 
suggesting that the transit is an eclipsing white dwarf- M 
dwarf binary, later confirmed by IMuirhead et al.l (|2013[ ) . 
The remaining sample of KOIs contains 157 planet can- 
didates orbiting 106 dwarfs. 

2.2. Kepler non-KOI sample 

IMann et all ()2012D show that > 90% of the bright 
{Kp < 14), and ~ 7% of the faint {Kp > 14), 
late K and M {Kp — J > 2.0) Kepler targets are gi- 
ant stars. JHK col ors are sometimes use d to iden- 
tify giant stars (e.g.. ILepine fc GaidosI l20Tll) . however, 

* http://exoplanetarchivc.ipac.caltech.edu/ 
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these colors are known to b e metal-sensitive (jLeggettl 
Il992t iMuirhead et al.|[2012b[ ). and the giant and dwarf 
branches overlap in J HK color space for K type stars 
(jBessell fc BrettI figSSt ). Instead, we screen out giants 
stars using their reduced proper motion, defined as: 



Hj = J + 51og/i + 5, 



(1) 



where /x is the total proper motion in arcseconds yr~^ 
and the J magnitude is taken from 2MASS. 

We computed proper motions for each Kepler target 
star using all available astro metry from the U SNO-Bl.O 
(JMonct ct al. 2003), 2MASS (iSkrutskie et alj f2006). and 
SDSS (Ahn ct al. 2012) catalogs, using the algorithm de- 
scribed in Kraus fc H illenb rand (2007) . We obtained the 
astrometry for each star from t he Vizier archive using 
the IDL routine queryvizier.pro ()Landsmanlll993D . and 
then combined the astrometry epochs from all surveys 
using a weighted least-squares fit. Our algorithm tested 
the goodness of each fit for each proper motion and re- 
jected all astrometry outliers at > Scr. Most of these 
outliers were found in the photographic survey data, not 
in 2MASS or SDSS, due to the heavy weight assigned to 
the modern CCD-based epochs. The resulting catalog 
has proper motion uncertainties a factor of ~2 smaller 
than those from USNO-B alone (~3 mas yr~^ versus 6-7 
mas yr^"'^). 

For bright [R < 12) stars, many of which are saturated 
in one or more of the above surveys, we adopted proper 
motions f rom the Third USNO CCD Astrograph Catalog 
(UCACa. lZacharias et al.][2010l) . UCAC3 extends to R = 
16, though the proper motion errors become quite large 
at i? > 13-14. The typical errors in the UCAC3 proper 
motions are ~l-3 mas yr~^ for stars as faint as i? ~ 12 
and ~6 mas yr~^ for those as faint as i? ~ 16. 

We use s tars with known luminosity class from 
iMann et al.l (|2012[ ) to test possible reduced proper mo- 
tion cuts. Figure [1] shows the reduced proper motions 
for the giant and dwarf samples, excluding those stars 
for which the errors in total proper motion are > 25 mas 
yr~^ and those that had potential contamination from a 
nearby star. A proper motion cut of Hj > 7.5 excludes 
only 1 dwarf (of 52) in the sample, and includes only 3 
giant stars (of 278). 

To establish the metallicity of Kepler target late K 
and M dwarfs, we randomly selected 100 stars observed 
by Kepler in Quarters 1 through 8 that have no detected 
planets (non-KOIs) for NIR spectroscopy with the cri- 
teria Kp — J > 1.85 and Hj > 7.5 (see Section [3] for a 
description of the observations) . Spectra of three of these 
100 contain st rong CO features at ^2.35 p,m indicativ e 
of giant stars (|Cushing et al.l 120051: IRavner et all 120090 . 
and were excluded from our analysis. 

3. OBSERVATIONS & REDUCTION 

Spectra of KOIs were obtained wi th the SuperNova 
Integral Field Spectrograph (SNIFS. Eantz et al.|[200l 
on the University of Hawaii 2.2m telescope atop Mauna 
Kea. SNIFS covers 3200 A to 9700 A at a resolution 
of 1000 < R < 1300. SNR was > 80 redward of 
6000 A for each ta rget. The SNIFS processing pipeline 
(jBacon et al.ll200"TI ) automatically performed basic data 
reduction. This included dark, bias, and fiat-field correc- 
tions, removing of bad pixels and cosmic rays, and sky 
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Figure 1. Distribution of reduced proper motions for spectro- 
scopically con firmed late K and M giants (black) and dwarfs (grey, 
hashed) from lMann et ahl ||2012I 1. Stars with questionable proper 
motions (error in proper motion > 25 mas yr~^) were excluded. 
We utilized a reduced proper motion cut of Hj > 7.5 (shown as 
a vertical dashed line) to remove interloping giant stars from our 
non-KOI sample. 

subtraction. The SNIFS pipeline used arcs taken at the 
same position as the target to wavelength- calib r ate th e 
data. Spectrophotometric standards from lOkd ()1990[) , 
taken over the course of each night, were used in conjunc- 
tion with a model of the atmosphere above Mauna Kea 
(But on et al.l l2013[ ) to correct for instrument response 
and atmospheric extinction. We shifted the spectra in 
wavelength to the rest frames of their emitting star by 
cross correlating them with a similar s pectral type tem- 
plate from the Sloan Digital Sky Survey (jBochanski et al.l 
[20f)l . 

We obtained near-IR spectra of t he 100 non-KOI tar- 
gets using the SpeX spectrograph (jRavner et al.l |2003() 
attached to the NASA Infrared Telescope Facility (IRTF) 
on Mauna Kea. SpeX observations were taken in the 
short cross-dispersed mode using the 0.3" slit, which 
yielded a resolution of i? ~ 2000 from 0.8 to 2.4/im. 
SNR in the H- and X-bands was typically > 80. To 
correct for telluric lines, we observed an AOV-type star 
within 30 minutes of time and 0.1 airmass of the tar- 
get observation. To remove effects from large telescope 
slews, we obtained fiat-field and argon lamp calibration 
sequences after each AOV star. Spectra were extracted 
and reduced using the SpeXTool package, which per- 
formed flat-field correction, wavelengt h calibration, and 
sky subtraction (jCushing et al.ll20M ). Telluric correc- 
tions were co mputed from each AOV star using the xtell- 
cor package (jVacca et al.ll2003t ). and then applied to the 
relevant target spectra. We then placed each spectrum in 
the star's rest frame by cross-correlating it to a spectrum 
of a template star (of a similar spectral type) from the 
IRTF spectral library (jCushing et al.ll2005l: IRavner et al.l 
12001 . 

4. DETERMINATION OF [FE/H] AND R, 

M13 provide empirical calibrations for calculating 
[Fe/H] from indices in visible, J-, H-, or K- band spec- 
tra of late K and M dwarfs. However, M13 have higher 
SNR observations (SNR ~ 150) than those obtained for 
the much fainter targets observed in this program. As 
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a result, when we applied the calibrations from Ml 3 on 
features blueward of 6000 A (NIR calibrations are less af- 
fected) the resulting errors were large (> 0.1 dex) from 
measurement noise alone. 

To mitigate SNR errors, we took the calibrator (wide 
binary) sample of M13 and repeated their process of 
defining a visible wavelength metallicity calibration. 
However, we restricted ourselves to indices redward of 
6000 A where the SNR of our observations is the highest. 
We then derived the following calibration: 

[Fe/H] =0.68F1 + 0.53i^2 - 0.32F3 

-l.OColorl - 0.26, (2) 

where Color 1 is a temperature sensitive index from 
iHawlev et all (|2002t ). and Fl, F2, and F3 correspond to 
the equivalent widths of features at 8191-8225 A, 8860- 
8880 A, and 9179-9199 A, respectivel}0. We used the 
pseudo-continuum regions defined in M13. Figure [2] 
shows the primary star metallicity as a function of the 
K/M dwarf metallicity derived from Eqn. [51 Equation [5] 
has a root mean square error (RMSE) of 0.08 dex, and 
an adjusted square of the multiple correlation coefficient 
(i?ap) of 0.82, indicating it is roughly as reliable as the 
corresponding calibration from M13 (RMSE=0.07 and 

^ap=0.84). 
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Figure 2. Metallicity of the primary (FGK) dwarf from M13 ver- 
sus metallicity for the late K or M dwarf companion derived using 
Eqn. [21 The dashed line indicates equality. The binary sample 
covers the full range of metallicities and spectral types in our KOI 
and non-KOI sample. 

Metallicities of the non-KOI sample were calculated 
as the weighted means of the J-, H-, and JC-band cal- 
ibrations outlined in M13. Weights were based on the 
measurement errors in each band added in quadrature 
with the errors from each calibration (e.g., the J-band 
calibration has a significantly higher RMSE than the H- 
and K-haiid calibrations). 

Using metallicities from two different sources (visible 
and NIR) engenders the risk of systematic differences. 
However, the empirical relations we utilized are cali- 
brated using an identical set of wide binaries. As a 
check , we selected a samp le of 55 late K to mid-M dwarfs 
from ILepine et al.l (|2013D that have both visible wave- 

^ Note: all wavelengths in this work are reported in vacuum 
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KOI Stellar Parameters 






KOI 


Kepler ID 


[Fe/H] 


<^[Fc/Hl 


Rq 


Rq 


227 


6185476 


-0.21 


0.08 


0.60 


0.03 


247 


11852982 


+0.11 


0.08 


0.53 


0.04 


248" 


5364071 


+0.13 


0.05 


0.57 


0.04 


249 


9390653 


-0.30 


0.08 


0.37 


0.06 


250" 


9757613 


-0.10 


0.08 


0.59 


0.03 


251" 


10489206 


-0.03 


0.08 


0.48 


0.05 


252 


11187837 


-0.08 


0.09 


0.50 


0.05 


253 


11752906 


+0.41 


0.07 


0.55 


0.04 


254 


5794240 


+0.27 


0.08 


0.49 


0.05 


255 


7021681 


-0.23 


0.09 


0.55 


0.04 


314" 


7603200 


-0.07 


0.08 


0.52 


0.04 


430 


10717241 


-0.12 


0.09 


0.64 


0.03 


463 


8845205 


-0.35 


0.08 


0.36 


0.06 


478 


10990886 


+0.20 


0.08 


0.48 


0.05 


503 


5340644 


-0.11 


0.09 


0.62 


0.03 



Note. — Table [T] is published in its entirety in 
the electronic edition of the Astrophysical Journal, 
and can be downloaded with the arXiv version of 
the manuscript. A portion is shown here for guidance 
regarding its form and content. 

'^ Multi-plane t cand idate system as listed in the 
IBatalha et"all 12013 ) catalog. 



length spectra from SNIFS and NIR spectra from SpeX. 
The scatter between [Fe/H] values derived using visi- 
ble and NIR spectra are consistent with combined mea- 
surement and calibration errors for these two techniques 
(a = 0.14 dex). More importantly, there is no significant 
offset between metallicities derived from the two meth- 
ods (median difference in metallicity is 0.01 ± 0.02 dex). 

We compa re our metallicit ie s of ou r KOI sample with 
values from [Muirhead et al.l (|2012af) in Fig. ^. Our 
metallicities are consistent with their's for stars with 
[Fe/H] > —0.3. However, for more metal-poor stars, 
metallicities from our analysis are systematically lower 
and the scatter between metallicity esti mates is higher. 
Note th at we report [Fe/H] values while iMuirhead et al.l 
(12012a) uses [M/H]. Increasing [a/Fe] with decreasing 
[Fe/H] may, in part, explain t his discrepancy. T he cali- 
l3rator (w ide binary) sample oflRoias-Avala et al.l ([2012), 
on which IMuirhead et al.l (|2012af) is based, has a paucity 
of stars with [Fe/H] < —0.5. The result is that their 
calibration assigns metallicities that are systematically 
too high for stars with [Fe/H] < —0.4. As can be seen 
in Fig. [2l our calibration performs well for dwarfs even 
with [Fe/H] < -0.5. 

In Fig. [5^) we c o mpare metallicities from 
[Dressing &: Charbonnead (|2013D with ours for 
all overlapping (KOI and non-KOI) targets. 
[Dressing fc Charbonnead ([2013[ ) [Fe/H] values are 

{■^ > 3), and show no 
[Dressing fc Charbonnead 



inconsistent with (reduced xl 

correl ation with our values. [I 

([2013() themselves note significant dis agreement with 

their m etallicities and those reported in [Muirhead et all 

(|2012a|) . highlighting the difficulties of estimating M 

dwarf metallicities from photometry and stellar models 

alone. 

We calculated s tellar radii usi ng the Tcff-i?* rela- 
tionship given in [Bovaiian et al.l ^2012). We deter- 
mined T^s for each KOI by fitting BT-SETTL models 
([Allard et al.,i2011,) to our visible wavelength spectra fol- 
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Table 2 

non-KOI Stellar Parameters 



Kepler ID 


[Fc/H] 


<^[Pc/Hl 


1721911 


+0.37 


0.07 


1996399 


+0.07 


0.06 


2010738 


+0.21 


0.12 


2850521 


-0.68 


0.11 


3233853 


-0.01 


0.05 


3342894 


+0.17 


0.09 


3533220 


-0.19 


0.06 


3935942 


+0.04 


0.05 


4543236 


-0.19 


0.06 


4543619 


-0.55 


0.12 


4553205 


+0.08 


0.05 


4682420 


-0.39 


0.12 


5000970 


-0.26 


0.07 


5165017 


-0.02 


0.06 


5252367 


-0.18 


0.13 



Note. — Table [2] is pub- 
lished in its entirety in the elec- 
tronic edition of the Astrophys- 
ical Journal, and can be down- 
loaded with the arXiv version 
of the manuscript. A portion is 
shown here for guidance regard- 
ing its form and content. 



lowing the technique outlined in iLepine et al.l ()2013D . 
except that we only included models with metallici- 
ties < 2cr different fro m thos e we d erived from Eqn. [21 
Stellar radii in iBatalha et al.l (120131) are b ased on tem- 
peratures in the KIC (iBrown et al.l 120 lit) and Yonsei- 
Yale isochrones (iDemarque et al.l I2004D . However, radii 
from IDemarque et al.l ( 20041) are known to be i naccu - 
rate for late K and M dwarfs (jBoyaiian et al.l I2012I ). 
and KIC temperatures ha ve been shown to b e too high 
for late K and M dw arfs (jMann et al.ll2012| ). Instead, 
iBovaiian et al.l (|2012D derive their Teg-R^ using empir- 
ical measurements of radii of K and M dwarfs from 
long-baseline interferometry. They obtain [Fe/H] val- 
ues for their targets from the literature, and find no 
discernible metallicity dependence in their T^g-R^, rela- 
tio n. This cont r adicts stellar evolutionary models such 
as iDotter et al.l ()2008[ ). which show a strong metallic- 
ity dependence for the Tes-R* relation fo r late K and M 
dwarfs. We chose to use the relation f r omlBovaiian et al.1 
( I201J), rather than t h e IDotter et all (I2008D models, be- 
cause IBovaiian et al.l (J2012D is based on empirical mea- 
surements rather than evolutionary models. 

Our stellar r a dii are 0.06i?o larger than those of 
iMuirhead et al.l ()2012aD . Their stellar radii for stars 
with Tcff< 3900 are consistent with ours at < 1.2cr. 
The discrepancy is larger (typically > 1.5cr) for warmer 
stars because their tem peratures are un derestimated: the 
II2O index utilized bv IMuirhead et al] (l2'012a) to calcu- 
late Te ff saturates at 3800-400(3 K (iRoias-Avala et al.l 
20121: iDressing fc Charbonneaul I2Q13D . However, 

Muirhead et al.l ()2012aD stellar radii are still < 3cr 
consistent even at warm temperatures. Our stel- 
lar r^du__aje__on__averageO .02^0 larger than those 
of IDressing fc Charbonneaul (|2013[ ) with no significant 
trend in Tgg. This offset is smaller than the typical errors 
from our own measurements (median (Tu^ ~ O.OARq). 
Resulting Tcs, -R*, and associated errors for KOIs in this 
program are listed in Table [Tj 
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Figure 3. Metallicities from IMuirhead et al.l Ij2012al ') (top) or 
from Dressing & Charbonncau (2013) (bottom) as a function of 
those derived from our own program. The dashed lines indicate 
equality. We have added artificial scat ter (~ 0.01 dex) to metal- 
licities from iPressing & Charbonncau 1)20131 ) for clarity. Typical 
errors are shown in th e bottom right of each plot. Note that 
IMuirhead et al.l Il2012al) use [M/H] instead of [Fe/H]. Our metal- 
licities are mostly within Icr (reduced x^ — t) of those from 
Muirh ead et al. (2012a), with the exception of those with [Fc/H] < 
—0.3 which are more discrepant. Our [Fe/H] values greatly differ 
(reduced x^ > 3) from those of lDressine fc Charbo nncau ( 20131) . 

We adopte d the Rvianet/R* values reported by 
IBatalha et all (|2013D . IDressing fc Charbonnea"ul (|2013l ) 
refit the Kepler light curves for M dwarf KOIs and 
find that there are significant proble ms with some 
of the transit parameters reported in IBatalha et al.l 
(|2013D . However, these problems are primarily in 
a/i?* (where a is the semi-major axis) and impact pa- 
rameter, whereas the niedian Rpianet/R* value from 
IDressing fc Charbonneaul (I2013D is only 3% smaller than 
those of IBatalha et al.l (I2013D. Further, the differences 
in R r,ianet/R* betw e en iDrcssing fc Charbonneaul (|2013f ) 
and IBatalha et al.l (|2Q13D are s mall compared to er - 
rors in Rpianet/R* reported by IBatalha et al.l (|2013f ) 
(c± 13%) and errors in stellar r adii (~ 7%). Moreover 
IDressing fc Charbonneaul (|2013| ) only refit transits of M 
dwarfs, while our sample includes many late K dwarfs. 

5. RESULTS 

We list [Fe/H] values and stellar radii for the KOI sam- 
ple in Table [1] and [Fe/H] values for the non-KOIs in 
Table [2] Fig. |4| compares the metallicity distributions 
of dwarfs with no detected transit (non-KOI sample), 
as well as metallicities for different planet-size samples 
(Earths, Neptunes, and Jupiters). We summarize the 
metallicities for each distribution in Table jS] and com- 
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Table 3 

Summary of Metallicity Distributions 



Planet- type 


Planet sizes 


N 


[Fe/H] 


KS" 


CMA'' 




R(B 




median 




a 


Jupiters 


Rp>8 


3 


+0.12 ±0.07 


2.0% 


3.2 


Neptunes 


2.0 <Rp <8 


40 


-0.10 ±0.04 


96.9% 


0.2 


Neptunes2 


2.5 <Rp <8 


17 


-0.12 ±0.06 


99.5% 


0.1 


Neptunes3 


3.0 <Rp <8 


6 


-0.03 ±0.11 


99.5% 


0.2 


Earths 


Rp<2 


63 


-0.15 ±0.03 


33.9% 


1.2 


Multis= 


all 


31 


-0.10 ±0.04 


39.4% 


0.8 


Non-KOI 


not detected 


97 


-0.10 ±0.03 







'^ Probability that the distribution is drawn from the same parent 
population as the non-KOI sample based on a Kolmogorov-Smirnov 
(KS) test. 

Difference between distribution and non-KOI distribution (in stan- 
dard deviations) as determined by the Cramer-von Mises- Anderson 
statistic (And erson 1962). 

'^ Systems wi t h mo re than one transiting planet detected in the 
IBatalha et al.l II2013I ) cataolog. 

pare with the non-KOI sample using the Kolmogorov- 
Smirnov (KS) and Cramer-von Mises- Anderson (CMA) 
tests. Systems with multiple detected planets are placed 
into size bins (Earth-, Neptune-, and Jupiter-size) ac- 
cording to the largest detected planet in the system. We 
also list the metallicity dis tribution of these mu lti-planet 
systems, according to the IBatalha et alj ()2013[ ) catalog, 
in Table O 

Metallicities of stars hosting Jupiter-sized planets are 
significantly higher than those with no detected transit, 
which is consistent with previous fi ndings based on radia l 
velocity surveys of M dwarfs (e.g., I Johnson et al.|[20Tol ). 
Because there are only 3 giant planets in the sample, the 
difference between the Jupiter and non-KOI sample is 
at the edge of statistical significance, with a difference 
between the median [Fe/H] values of 0.22 ±0.073 (3. Oct). 

The distribution of metallicities for dwarfs hosting 
Neptune-sized planets is consistent with the non-KOI 
sample for all metrics. We rule ou t the .2 dex offset 
reported bv lSchlaufman fc LaughliiJ (|2011D between the 
metallicity of stars hosting no n- Jovian (Rp < SRq) plan- 
ets and the non-KOI sample at > 6cr, and at > ia if we 
consider just Neptune-sized planets (2.0 < Rp < SR®). 

The distribution for the Earth-sized hosts is slightly 
more metal poor (by 0.05 dex) than the non-KOI sam- 
ple, although the offset is not significant {a = 0.04 dex). 
Our detected offset is consistent with predictions from 
iGaidos fc MannI (|2013[) . who use iDotter et al.l ()2008f ) 
models and Kepler target stars to show that M dwarfs 
hosting small planets have [Fe/H] values ~ 0.02 dex lower 
than those without planets, because for a given g — r 
color, metal-poor K and M dwarfs will have smaller radii 
than metal-rich dwarfs (maki ng small planets easie r to 
detect). If the relations from iBovaiian et alJ (|2012l ) are 
correct (i.e., metallicity is negligible factor in i?* for a 
given Teff) th en this detection bias is smaller, but still 
present, since IBovaiian et all (|2012[ ) find that the color- 
i?* and color- Tcff relations for K and M dwarfs have a 
significant metallicity dependence. 

6. DISCUSSION 

In this paper we present our comparison between 
metallicities of late K and M Kepler target stars and 
planet candidate hosts. We used techniques of (or 



modified techniques of) iMann et al.l (|2013D to calculate 
[Fe/H] . We then investigated correlations between stellar 
metallicity and the presence, multiplicity, and size of any 
detected planets. We draw four main conclusions: 

• The metallicity distribution of late K and M Kepler 
targets is indistinguishable from that of the solar 
neighborhood. 

• Late K and M Kepler dwarfs hosting giant planets 
are more metal-rich than those without detected 
planets. 

• Late K and M hosts where the largest detected 
planet is Earth- or Neptune-sized have metallici- 
ties consistent with those dwarfs with no detected 
transit. 

• Late K and M dwarfs hosting multiple detected 
planets are not significantly more metal-rich or 
metal-poor than those with no detected transit. 

An important effect is the presence of non-detected 
planets in the control (non-KOI) sample, which dilutes 
any metallicity offset between the two samples. Suppose 
metallicity is a bimodal distribution with stars harboring 
planets having metallicity greater by A[fo/h] than those 
without planets. The observed metallicity offset (O) be- 
tween the KOI sample and the non-KOI sample is: 



= A 



[Fc/H] 



(1-/) 

{i-fty 



(3) 



where, / is the fraction of stars with planets, and t is 
the probability of detecting the planet (e.g., the geo- 
metric transit probability). For transiting planets, ft is 
small, and the denominator ~ 1. In the case that / ap- 
proaches 1, O ~ 0, because the non-KOI sample is com- 
pletely diluted with undetected planets. This may be the 
case when considering Earth- to super-Earth -sized plan- 
ets an d all orbital periods around M dwarfs () Swift et alJ 
|2013[) . In the case of giant plan ets, which are relatively 
rare (2% for period <85 davs: jPVessin et ahl ^20131. di- 
lution is negligible and O ~ Apo/H]- For Neptunes, 
iFressin et al.l (|2013f) find that ~25% of stars harbor a 
Neptune (2i?0 < Rp < 8 Rq) and periods <85 days 
(note that only 1 of the planets in our sample has a pe- 



0.75A 



[Fo/H]' 



which has 



riod 3>85 days). In this case O 
little effect on our conclusions 

We examined how our results change as a function of 
how we define Earth-, Neptune-, and Jupiter-sized plan- 
ets. There are no KOIs in our sample with radii between 
6i?0 and 9i?0 , so our choice of a Neptune- Jupiter bound- 
ary is unimportant. We investigated the effect of chang- 
ing the Earth- Neptune boundary by considering two sub- 
samples; Neptunes2, defined as 2.5 R^ < Rp < 8 i?©, and 
Neptunes3, defined as 3.0 i?© < Rp < 8 /J®. All three 
Neptune samples are consistent with each other and the 
non-KOI sample at Icr, demonstrating that our results 
are not sensitive to how we define Earth- and Neptune- 
si zed planets. 

IFressin et al.l ()2013l ) show that the Kepler planet search 
algorithm is missing planets that should have been de- 
tected based on their SNR. However, this only alters our 
results if the metallicity distribution of missing planets 
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Figure 4. Metallicity distribution of Kepler late K and M targets with no detected transit (black), Earth-sized KOIs {Rp < 2i?0, 
magcnta-grey dashed), Neptune sized KOIs {2Rq < Rp < SiJg, blue), and Jupiter-sized KOIs {8Rq < Rp, red, hashed). Bin sizes and 
locations are identical for all distributions, but the histograms are slightly offset for clarity. Arrows above the plot indicate the median of 
each distribution. Note that the black and blue arrows are nearly overlapping. For multi-planet systems we use the radius of the largest 
detected planet in the system. 



is significantly different than that of the detected plan- 
ets. To check the effect of incompleteness, we consid- 
ered a subsample with transit detections of SN R > 16 
as rep orted by Batalha et al. (2013), where Fressi n et al.l 
()2013() suggest the detection efficiency of the Kepler 
pipeline is ~ 100%. In this case, the sample shrinks 
from 157 planet candidates around 106 dwarfs to 132 
candidates around 93 dwarfs. Most of the candidates re- 
moved by this cut are Rp < 2Rq. We reran all analyses 
on this subsample and find that none of our conclusions 
are changed. 

Our results di s agree with those of 
iSchlaufman fc LaughlinI (|2011[ ). who claimed that late 
K and M Kepler stars hosting small planets are more 
metal-rich than non-hosting late K and M stars. In place 
of sp ectroscopic metallicities Schlaufman & Laughlin 
()2011|) use g — r versus J ~ H colors, which been shown 
to be positively correlated with metallici ty for F, G, 
and early-K dwarfs ('e. g.. lAn et al.l I2009D . However, 
iWest et all (|2004f ) and iBochanski et al.l (|2013l ) see a 
negative correlation between g — r color and metallicity 
for M dwarfs. The stars in our sample arc in the 
transition region (mostly K5-M2), where the PHOENIX 
stellar atmosphere models predict little or no trend of 
g — r color with metallicity (jLepine et al.ll201^ . 

We used our sample (both KOIs and non-KOIs) to in- 
vestigate how g — r correlates with metallicity for our 
range of spectral types. In Fig. [5^ we show the g — r 
colors for the metal-poor ([Fe/H] < —0.1) and metal- 
rich ([Fe/H]> —0.1) samples for three different J — H 
bins centered at J — iJ = 0.575, 0.625, and 0.675. 



[Fe/H]= —0.1 was selected to divide the samples because 
this is the median me tallicity of our sample. The la te- 
type dwarf bin used bv ISchlaufman fc LaughlinI (|2011| ) is 
centered at J — H — 0.62. The distribution oi g — r colors 
of the two metallicity samples are consistent at < 2ct in 
each of the three J — H bins. We use a slightly different 
parsing in Fig. [Sb, where we show [Fe/H] versus g — r 
color for two J — H bins. In each bin a least-squares lin- 
ear fit to the data yield slopes that are not significantly 
different from 0. The fit yields coefficients of determi- 
nation (R2) of 0.01 and 0.07 for the 0.65 > J - 7J and 
0.65 < J — H bins, respectively. An F-test comparing 
[Fe/H] -[Fe/H] to that of [Fe/H]- [Fe/H] /,t gives respec- 
tive differences in the variances of only 0.50 and 0.48cr. 
These results strongly suggest that g — r versus J — H is 
not a good predictor of [Fe/H] for late-type dwarfs. 

Our results have some important theoretical impli- 
cations. Theoretical studies have suggested there is a 
minimum me tallicity for a protop l anetary disk to form 
planets fe.g.. iGonzalez et al.ll200lt iJohnson fc Lill2012| ). 
Assuming the metallicity of the protoplanetary disk 
matches that of the star later in its evolution, this sug- 
gests thB;t_Jow;^nBtalhcity stars should not harbor plan- 
ets. IJohnson fc Lil ()2012D estimate that the minimum 
metallicity to form a planet around a solar-type star is 
[Fe/H]niin ~ — 1.5-1- log(a), where a is the semi-major axis 
in AU. For the planets in our paper (a < 0.1) [Fe/H]min 
is approximately —2.5. Because M dwarfs have less mas- 
sive, we expect that [Fe/H]min will be higher for these 
stars. Our results show that Earth and Neptune-sized 
planets are able to form around stars with metallicities 
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Figure 5. Left: Median g — r colors of the metal-rich (blue) and metal-poor (red) samples of Kepler dwarfs for three J — H color bins 
centered at 0.575, 0.625, and 0.675. Bin sizes and locations are identical for all distributions, but are slightly offset for clarity. The scatter 
in g — r colors for a given inetallicity range and J — H bin are determined through bootstrap resampling. For all three J — H bins, the two 
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both bins, the best fit lines are consistent with a slope of 0, and yield coefficients of determination (R^) of 0.01 and 0.07, respectively. An 
F-tcst detects no significant improvement from the regression for either bin. Thus no significant correlation between g — r versus J — H 
and [Fe/H] is found. 



as low a s [Fe/Hjc::; —0.6, simila r to what is seen for FGK 
dwarfs ()Bucliliave et al.ll2012[ ). But it is likely that we 
are not probing sufficiently metal-poor dwarfs to detect 
the proposed planet formation threshold. 

Our results also indicate that, for small planets, mul- 
tiplicity is not correlated with metallicity. Interestingly, 
two of the four KOIs with [Fe/H] < -0.5 have > 3 de- 
tected planets (KOI-961 has 3 and KOI-812 has 4), sug- 
gesting the accretion process must be efficient in collect- 
ing solids from the disk. A minimum mass solar neb- 
ula contains about 64M0 of rock/m etal/ices. Assum- 
ing t hat the disk mass is ^ 0.1 A/* ([Williams fc Ciezal 
1201 If) , and that the amount of metals in a disk scales 
with its mass and the stellar metallicity, a disk around 
an early M dwarf with [Fe/H] = —0.5 contains ^ lOM^ 
of solids. KOI-812 (as an example) contains 4 detected 
planets with radii from 1.3 to 2.4R(^. Most of these are 
likely rocky, or are c omposed of rocky c ores with a thin 
hydrogen envelopes (jGaidos et al.ll201^ . If we assume a 
mass-radius relationship of Mp ~ Rj^, with Mp and Rp 
in Earth units, then the total mass in KOI-812's planets 
is 15.5 ± 3.1M0. Although some of this mass is hydro- 
gen (and thus not affected by the amount of solids in the 
disk), this analysis does not consider undetected planets 
at higher semi-major axes. Thus our results suggest that 
the progenitors of these planets must have been very ef- 
ficient in accreting most of the available solid material 
from the disk. 

The core accretion scenario of giant planet forma- 

tion requires the formation of a. ^ 5 — lOM^ core 

()Pollacket al.l [19961 : iHubickvi etaD [20051 ) in the ~ 2 - 
6 Myr time scale on which disks are obs erved to dissipate 
([Haisch et al. 2001: lEvans et al.ll2009t) . The scarcity of 
giant planets around M dwarfs (especially metal-poor M 
dwarfs) means that either giant planets cores do not form 
around these stars, or that they do not form in time. 
The existence of objects that are likely 5 — lOM^ (or 
have rocky cores of this size), even around metal-poor M 
dwarfs, suggests that the latter is the more viable expla- 
nation. 



We thank the anonymous reviewer for helping to make 
the manuscript significantly better. This work was 
supported by NSF grant AST-0908419, NASA grants 
NNX10AI90G and NNX11AC33G (Origins of Solar Sys- 
tems) to EG. 

This paper includes data collected by the Kepler mis- 
sion. Funding for the Kepler mission is provided by the 
NASA Science Mission Directorate. SNIFS on the UH 
2.2-m telescope is part of the Nearby Supernova Factory 
project, a scientific collaboration among the Centre de 
Recherche Astronomique de Lyon, Institut de Physique 
Nuclaire de Lyon, Laboratoire de Physique Nuclaire et 
des Hautes Energies, Lawrence Berkeley National Labo- 
ratory, Yale University, University of Bonn, Max Planck 
Institute for Astrophysics, Tsinghua Center for Astro- 
physics, and the Centre de Physique des Particules de 
Marseille. Based on data from the Infrared Telescope Fa- 
cility, which is operated by the University of Hawaii un- 
der Cooperative Agreement no. NNX-08AE38A with the 
National Aeronautics and Space Administration, Science 
Mission Directorate, Planetary Astronomy Program. 

Facilities: IRTF, UH:2.2m, Kepler 

REFERENCES 



Adams, F. C, Bodenheimer, P., & Laughhn, G. 2005, 
Astronomische Nachrichten, 326, 913 

Ahn, C. P., et al. 2012, ApJS, 203, 21 

AUard, F., Homeier, D., & Freytag, B. 2011, in Astronomical 
Society of the Pacific Conference Series, Vol. 448, 16th 
Cambridge Workshop on Cool Stars, Stellar Systems, and the 
Sun, ed. C. Johns-KruU, M. K. Browning, & A. A. West, 91 

An, D., et al. 2009, ApJ, 707, L64 

Anderson, T. W. 1962, Annual Math Statistics, 33, 1148 

Bacon, R., et al. 2001, MNRAS, 326, 23 

Batalha, N. M., et al. 2010, ApJ, 713, L109 

— . 2013, ApJS, 204, 24 

Bessell, M. S., & Brett, J. M. 1988, PASP, 100, 1134 

Bochanski, J. J., Savcheva, A., West, A. A., & Hawley, S. L. 
2013, AJ, 145, 40 

Bochanski, J. J., West, A. A., Hawley, S. L., & Covey, K. R. 
2007, AJ, 133, 531 

Borucki, W. J., et al. 2010, Science, 327, 977 

Boyajian, T. S., et al. 2012, ApJ, 757, 112 



Metallicities of Kepler M dwarfs 



Brown, T. M., Latham, D. W., Everett, M. E., & Esquerdo, G. A. 

2011, AJ, 142, 112 
Buchhave, L. A., et al. 2012, Nature, 486, 375 
Buton, C, et al. 2013, A&A, 549, A8 

Casagrande, L., Flynn, C, & Bessell, M. 2008, MNRAS, 389, 585 
Gushing, M. G., Rayner, J. T., & Vacca, W. D. 2005, ApJ, 623, 

1115 
Gushing, M. G., Vacca, W. D., & Rayner, J. T. 2004, PASP, 116, 

362 
Demarque, P., Woo, J.-H., Kim, Y.-G., & Yi, S. K. 2004, ApJS, 

155, 667 
Dhital, S., West, A. A., Stassun, K. G., Bochanski, J. J., Masscy, 

A. P., & Bastien, F. A. 2012, AJ, 143, 67 
Dotter, A., Chaboyer, B., Jevremovic, D., Kostov, V., Baron, E., 

& Ferguson, J. W. 2008, ApJS, 178, 89 
Dressing, C. D., & Charbonneau, D. 2013, ApJ, 767, 95 
Evans, N. J., et al. 2009, ApJS, 181, 321 
Fabrycky, D. C, et al. 2012, ArXiv e-prints 
Fischer, D. A., & Valenti, J. 2005, ApJ, 622, 1102 
Fressin, F., et al. 2013, ApJ, 766, 81 
Gaidos, E., Fischer, D. A., Mann, A. W., & Lepinc, S. 2012, ApJ, 

746, 36 
Gaidos, E., & Mann, A. W. 2013, ApJ, 762, 41 
Gonzalez, G. 1997, MNRAS, 285, 403 

Gonzalez, G., Brownlee, D., & Ward, P. 2001, Icarus, 152, 185 
Haisch, Jr., K. E., Lada, E. A., & Lada, G. J. 2001, ApJ, 553, 

L153 
Hartmann, L., Calvet, N., Gullbring, E., & D'Alessio, P. 1998, 

ApJ, 495, 385 
Hawley, S. L., et al. 2002, AJ, 123, 3409 
Howard, A. W., et al. 2012, ApJS, 201, 15 
Hubickyj, O., Bodenheimer, P., & Lissauer, J. J. 2005, Icarus, 

179, 415 
Johnson, J. A., AUer, K. M., Howard, A. W., & Crepp, J. R. 

2010, PASP, 122, 905 
Johnson, J. A., & Apps, K. 2009, ApJ, 699, 933 
Johnson, J. L., & Li, H. 2012, ApJ, 751, 81 
Kokubo, E., Kominami, J., & Ida, S. 2006, ApJ, 642, 1131 
Kraus, A. L., & Hillenbrand, L. A. 2007, AJ, 134, 2340 
Landsman, W. B. 1993, in Astronomical Society of the Pacific 

Gonference Series, Vol. 52, Astronomical Data Analysis 

Software and Systems II, ed. R. J. Hanisch, R. J. V. 

Brissenden, & J. Barnes, 246 
Lantz, B., et al. 2004, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 5249, 

Society of Photo-Optical Instrumentation Engineers (SPIE) 

Gonference Series, ed. L. Mazuray, P. J. Rogers, & 

R. Wartmann, 146-155 



Laughlin, G., Bodenheimer, P., & Adams, F. G. 2004, ApJ, 612, 

L73 
Leggett, S. K. 1992, ApJS, 82, 351 
Lcpine, S., & Gaidos, E. 2011, AJ, 142, 138 
Lcpine, S., Hilton, E. J., Mann, A. W., Wilde, M., Rojas-Ayala, 

B., Cruz, K. L., & Gaidos, E. 2013, AJ, 145, 102 
Lepine, S., Rich, R. M., & Shara, M. M. 2007, ApJ, 669, 1235 
Mann, A. W., Brewer, J. M., Gaidos, E., Lepine, S., & Hilton, 

E. J. 2013, AJ, 145, 52 
Mann, A. W., Gaidos, E., Lcpine, S., & Hilton, E. J. 2012, ApJ, 

753, 90 
Monet, D. G., et al. 2003, AJ, 125, 984 
Muirhead, P. S., Hamren, K., Schlawin, E., Rojas-Ayala, B., 

Govey, K. R., & Lloyd, J. P. 2012a, ApJ, 750, L37 
Muirhead, P. S., et al. 2012b, ApJ, 747, 144 
— . 2013, ApJ, 767, 111 
Oke, J. B. 1990, AJ, 99, 1621 
Pollack, J. B., Hubickyj, O., Bodenheimer, P., Lissauer, J. J., 

Podolak, M., & Greenzweig, Y. 1996, Icarus, 124, 62 
Rayner, J. T., Gushing, M. G., & Vacca, W. D. 2009, ApJS, 185, 

289 
Rayner, J. T., Toomey, D. W., Onaka, P. M., Denault, A. J., 

Stahlberger, W. E., Vacca, W. D., Gushing, M. G., &; Wang, S. 

2003, PASP, 115, 362 
Rojas-Ayala, B., Covey, K. R., Muirhead, P. S., & Lloyd, J. P. 

2010, ApJ, 720, L113 
— . 2012, ApJ, 748, 93 

Santos, N. C, Israehan, G., & Mayor, M. 2001, A&A, 373, 1019 
Schlaufman, K. C, & Laughlin, G. 2010, A&A, 519, A105 
— . 2011, ApJ, 738, 177 

Scholz, A., Jayawardhana, R., & Wood, K. 2006, ApJ, 645, 1498 
Skrutskie, M. F., et al. 2006, AJ, 131, 1163 
Sousa, S. G., et al. 2008, A&A, 487, 373 
Swift, J. J., Johnson, J. A., Morton, T. D., Crepp, J. R., Montet, 

B. T., Fabrycky, D. G., & Muirhead, P. S. 2013, ApJ, 764, 105 
Terrien, R. C, Mahadevan, S., Bender, C. F., Deshpande, R., 

Ramsey, L. W., & Bochanski, J. J. 2012, ApJ, 747, L38 
Vacca, W. D., Gushing, M. C, & Rayner, J. T. 2003, PASP, 115, 

389 
West, A. A., et al. 2004, AJ, 128, 426 
Williams, J. P., & Gieza, L. A. 2011, ARA&A, 49, 67 
Wolfgang, A., & Laughlin, G. 2012, ApJ, 750, 148 
Woolf, V. M., Lepine, S., & Wallerstein, G. 2009, PASP, 121, 117 
Zacharias, N., et al. 2010, AJ, 139, 2184 



